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Abstract 

Effect of uniaxial single-ion anisotropy upon magnetic properties of a mixed spin- 
1/2 and spin-S" (S > 1) Ising model on a bathroom tile (4-8) lattice is examined 
within the framework of an exact star-triangle mapping transformation. The par- 
ticular attention is focused on the phase diagrams established for several values of 
the quantum spin number S. It is shown that the mixed-spin bathroom tile lattice 
exhibits very similar phase boundaries as the mixed-spin honeycomb lattice whose 
critical points are merely slightly enhanced with respect to the former ones. The 
influence of uniaxial single-ion anisotropy upon the total magnetization vs. temper- 
ature dependence is particularly investigated as well. 
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1 Introduction 



Over the last six decades, much effort has been devoted to determine a critical- 
ity and other statistical properties of various lattice-statistical models, which 
would enable a deeper understanding of order-disorder phenomena in solids. 
The planar Ising model takes a prominent place in the equilibrium statis- 
tical mechanics from this point of view as it represents a rare example of 
exactly solvable lattice-statistical model since Onsager's pioneering work [1]. 
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It is worthy to notice, however, that exactly soluble planar Ising models [2] 
are of immense interest because of providing a convincing evidence for many 
controversial results predicted in the phase transition theory and moreover, 
they also provide a good testing ground for approximative theories. Neverthe- 
less, it should be mentioned that a precise treatment of planar Ising models 
is usually connected with the usage of sophisticated mathematical methods, 
which consequently lead to considerable difficulties when applying them to 
more complicated but simultaneously more realistic models. 

The Ising systems consisting of mixed spins of different magnitudes, so-called 
mixed-spin Ising models, being among the most interesting extensions of the 
standard spin-1/2 Ising model which are currently of great research inter- 
est. These models have recently enjoyed a considerable attention chiefly due 
to much richer critical behavior they display compared with their single-spin 
counterparts. Magnetic properties of the mixed-spin Ising models can essen- 
tially be modified, for instance, by a presence of single-ion anisotropy acting 
on the S > 1 spins. Indeed, this effect may potentially cause a tricritical phe- 
nomenon - a change of the character of phase transition from a second-order 
to a first-order one. Another aspect which started to attract an appreciable 
interest towards the mixed-spin Ising models can be related to a theoretical 
modeling of magnetic structures suitable for describing a ferrimagnetism of 
certain class of insulating magnetic materials. In this respect, the ferrimag- 
netic mixed-spin Ising models are very interesting also from the experimental 
point of view, first of all in connection with many possible technological ap- 
plications of ferrimagnets to practice (e.g. thermomagnetic recording). 

Despite the intensive studies, there are only few examples of exactly solvable 
mixed-spin Ising models yet. Using the generalized form of decoration-iteration 
and star-triangle transformations, the mixed spin-1/2 and spin-S* (S > 1) 
Ising models on the honeycomb, diced and decorated honeycomb lattices were 
exactly been treated by Fisher long ago [3]. Notice that these mapping trans- 
formations were later on further generalized in order to account also for the 
single-ion anisotropy effect. It is worth mentioning that in this way general- 
ized mapping transformations were recently employed to obtain exact results 
of the mixed-spin Ising models on the honeycomb lattice [4] as well as on some 
decorated planar lattices [5]. To the best of our knowledge, these are the only 
mixed-spin planar Ising models with generally known exact solutions except 
several mixed-spin Ising models on the Bethe (Cayley tree) lattices studied 
within a discrete non-linear map [6] and an approach based on exact recursion 
equations [7]. Among the remarkable models for which a precise solution is 
restricted to a certain subspace of interaction parameters only, one should fur- 
ther mention the mixed-spin Ising model on the Union- Jack (centered square) 
lattice treated within the mapping onto a symmetric eight- vertex model [8]. 

The purpose of this work is to provide an exact formulation of the mixed 
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spin-1/2 and spin-.S (S > 1) Ising model on the bathroom tile (4-8) lattice 
and to establish accurate phase diagrams of this system for several values of 
the quantum spin number S. Exact results for the system under consideration 
are obtained by applying the generalized star-triangle mapping transformation 
constituting an exact correspondence with an effective spin-1/2 Ising model 
on the Shastry-Sutherland (orthogonal-dimer) lattice (further abbreviated as 
SSL), which has recently been solved by the present author elsewhere [9]. 
The obtained phase diagrams will also be compared with those of another 
exactly soluble three-coordinated planar Ising model, namely, on the mixed- 
spin honeycomb lattice [4]. Next, the influence of uniaxial single-ion anisotropy 
upon the temperature dependence of the total magnetization will be clarified. 

The outline of this paper is as follows. In Section 2, a detailed formulation of 
the considered model is presented and subsequently, exact expressions for the 
magnetization and critical temperatures are derived within the exact mapping 
onto the effective spin-1/2 Ising model on the anisotropic SSL. The most 
interesting numerical results obtained for several values of the quantum spin 
number S are presented and particularly discussed in Section 3. Finally, some 
concluding remarks are given in Section 4. 



2 Model and its exact solution 

Let us begin by considering a mixed spin-1/2 and spin-S* (S > 1) Ising model 
on the bathroom tile (4-8) lattice schematically illustrated on the left-hand 
side of Fig. 1. The mixed-spin 4-8 lattice consists of two interpenetrating 
inequivalent sub-lattices A and B, which are formed by the spin-1/2 (empty 
circles) and spin-S* (solid circles) atoms, respectively. The Ising Hamiltonian 
defined upon the underlying 4-8 lattice reads: 

3N N 

n^ 8 = -j £ sfsf-D^sf) 2 , (i) 

(i,j)cC i=l 

where Sf = ±1/2 and Sf = —S,—S + 1,...,S are the Ising spins located 
at ith lattice point, their superscript denotes the sub-lattice to which they 
belong, J being the exchange interaction between nearest-neighboring A — B 
spin pairs and the parameter D measures a strength of the uniaxial single-ion 
anisotropy acting on the spins of sub-lattice B. Furthermore, N denotes the 
total number of sites at each sub-lattice, the first summation is carried out 
over all nearest-neighboring A — B spin pairs and the second summation runs 
over all sites of sub- lattice B. 

It is of principal importance that the 4-8 lattice belongs to 'loose-packed' 
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lattices, i.e. it can be divided into two interpenetrating sub-lattices A and B 
in such a way that all nearest neighbors of A-sites belong to the sub-lattice 
B and vice versa. Owing to this fact, a summation over the spins of one sub- 
lattice (say sub-lattice B) can be performed before summing over the spins 
of another sub-lattice (sub-lattice A) when developing the partition function 
of the system under investigation. In the consequence of that, the partition 
function of the mixed-spin 4-8 lattice can be rewritten in this useful form: 

N +S 

Z ^= Ell E eMPD(S*) 2 }eM{3JS?(Sf + S£ + Sf)], (2) 

{S A } i=l Sf =-S 

where (3 = l/lksT), k-Q being Boltzmann's constant, T stands for the absolute 
temperature, the product is taken over all -B-sites and the symbol J2{s A } means 
a summation over all possible spin configurations of the sub-lattice A. It is 
quite obvious that this form justifies a performance of the familiar star-triangle 
mapping transformation [3]: 

£ e W [mS*) 2 ]eMf3JS?(Sf + S£ + Sf)] 
sf=-s 

= Cexp[/3J inter (SX + S^Sf + SfSf)}, (3) 

which establishes an exact correspondence between the mixed-spin Ising model 
on the 4-8 lattice and its equivalent spin-1/2 Ising model on the SSL when 
applying (3) to all -B-type sites (Fig. 1). Notice that the unknown mapping 
parameters C and Ji n t er are unambiguously given directly by the transforma- 
tion (3), which must hold for all available configurations of Sf, S^, spins. 
According to this, one readily finds: 

C = vl j Wl j \ p.J intCT = hy(Vi/V 2 ), (4) 

with the expressions V\ and V 2 defined as: 

s 

V 1= ]T exp(/? J Dn 2 )cosh(3/3Jn/2), (5) 

n=-S 
S 

V 2 = exp( y [3Dn 2 )cosh( y [3Jn/2). (6) 

n=-S 

At this stage, let us substitute the mapping (3) into the formula (2) gained 
for the partition function of mixed-spin 4-8 lattice. After straightforward re- 
arrangement, one easily obtains an exact relationship between the partition 
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function Z^s of the mixed-spin 4-8 lattice and its corresponding partition 
function Z SSL of the spin- 1/2 Ising model on the anisotropic SSL: 



Z 4 _ 8 (P, J, D) — C N Z S Sl(P, ^intra, Winter), 



(7) 



definitely given by the effective Hamiltonian: 




nj j \ ^ qA qA j \ ~* qA qA 

riSSL — — ^intra ~ Winter ' 




(ij)CV (ij)cC 



Above, the first (second) summation accounts nearest-neighbor intra-dimer 
(inter-dimer) interactions of the SSL that are displayed on the right-hand side 
of Fig. 1 as solid (dashed) lines and simultaneously, the condition Ji ntr a = 
2J inter must be satisfied. With regard to this, the isotropic mixed-spin 4-8 
lattice is mapped onto the anisotropic spin-1/2 Ising model on the SSL recently 
solved by the present author elsewhere [9]. The intra-dimer coupling of the SSL 
is twice as large as the inter-dimer one because the former interaction is being 
counted twice when performing the star-triangle transformation (3) at both 
opposite 5-sites of the same square plaquette of the original 4-8 lattice, while 
the latter one is being counted just once. 

The mapping relation (7) between the partition functions Z^ 8 and Zssl rep- 
resents a central result of our calculation, in fact, it can be utilized for estab- 
lishing similar mapping relations also for other important quantities such as 
Gibbs free energy, internal energy, magnetization, correlation functions, spe- 
cific heat, etc. When combining the formula (7) with commonly used mapping 
theorems [10], one easily proves a validity of following relations: 



where the symbols (...)4_8 and (...)ssl denote standard canonical average per- 
formed on the Ising lattices defined by means of the Hamiltonian (1) and (8), 
respectively, and the function F(x) labels expression: 



m A = (Sf) 4 -8 = (Sf)ssL = m S sL] 

i / oA oA oA\ _ / oA oA oA\ j. 

l A = Wj *Jk b l /4-8 — Wj /SSL = tSSL, 

m B = (Sf ) 4 _ 8 = 3m A [F(3) + F(l)]/2 + 2t A [F(3) - 3F(1)]; 



(9) 
(10) 

(11) 



s 



exp((3Dn 2 )nsmh(pjnx/2) 



F(x) 



n 



.=-S 



(12) 



s 



exp((3Dn 2 ) cosh((3Jnx/2) 



n=-S 
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It should be pointed out that the set of Eqs. (9)-(ll) enables straightforward 
calculation of both sub-lattice magnetization and of the mixed-spin 
Ising model on the 4-8 lattice. According to Eq. (9), the sub-lattice magne- 
tization rriA directly equals to the corresponding magnetization m SS L of the 
SSL being unambiguously given by the mapping relation (4) and the con- 
straint Jmtra = 2Ji ntcr - On the other hand, the triplet order parameter tssL is 
required for a computation of the sub-lattice magnetization mg in addition 
to m S sL- Although a precise evaluation of the triplet correlation t$sL between 
three spins located within the triangular plaquettes of the SSL has not been 
presented in Ref. [9], it can be accomplished following the same procedure as 
presented there for mssL- Subsequently, one readily proves that the triplet 
correlation tssL can be expressed in terms of the magnetization tjissl only. 
In view of these facts, an exact solution for both sub-lattice magnetization of 
the mixed-spin Ising model on the 4-8 lattice is thus formally completed. 

Finally, an analytical condition determining exact critical points of the mixed- 
spin 4-8 lattice will be derived. It can easily be understood that both sub- 
lattice magnetization and mg vanish when the magnetization mssL of the 
corresponding SSL tends to zero. In order to provide an accurate condition for 
the criticality of mixed-spin 4-8 lattice it is therefore sufficient to ascertain an 
exact critical point of the anisotropic SSL fulfilling the claim Jmtra = 2J inter . 
After straightforward but a little bit tedious modify based on the procedure 
presented in Ref. [9], one finds an exact critical point of the anisotropic SSL 
with this special ratio of inter- and intra-dimer coupling constants being: 



„ T Jintra 01 ^ 8 + 5^ + ^10 + 8^2 + 2^5 + 4y/2 \ 

Pc Jintra = T 7p~ = 2 hi . (13) 

ftBJc ' 2 + ^+^10 + 8^ 



Regarding this, a direct substitution of the critical point (13) into the relevant 
mapping relation (4) immediately yields the critical condition determining 
phase boundaries of the mixed-spin Ising model on the 4-8 lattice given by: 



s 

exp(P c Dn 2 ) cosh(3/3 c Jn/2) 

n=-S _ 

S ~ 

^ exp((3 c Dn 2 ) cosh(/3 c Jn/2) 

n=-S 

8 + 5^+^10 + 8^ + 2^5 + 4^2 



2 + ^+^0 + 8^2 



6 



3 Results and Discussion 



Let us step forward to a discussion of the most interesting numerical results. 
Although all results presented in the preceding section hold for the ferromag- 
netic (J > 0) as well as ferrimagnetic (J < 0) version of the model under 
investigation, in what follows we shall restrict ourselves to an analysis of the 
ferrimagnetic model only. It appears worthwhile to mention, nevertheless, that 
phase diagrams displayed below for the ferrimagnetic model are maintained 
without any changes also for the ferromagnetic model as a result of an invari- 
ance of the critical condition (14) with respect to the transformation J — > — J. 

The uniaxial single-ion anisotropy basically modifies the magnetic behavior of 
the spin-S* atoms that constitute the sub-lattice B of mixed-spin 4-8 lattice. 
To illustrate the case, our attention was firstly focused on the ground-state 
behavior explored for arbitrary quantum spin number S. When the strength 
of the easy-plane single-ion anisotropy (D < 0) reinforces, the system exhibits 
successive first-order phase transitions associated with lowering the spin mo- 
ment S B from its highest possible value. The cascade of the uniaxial single-ion 
anisotropies at which these first-order transitions occur is given by: 



rr = — • < 15 > 

\J\ 4n-2' K J 

where n — 3/2, 5/2, S for half-integer spin-S* atoms, while n = 1, 2, S for 
integer spin-S* atoms. The most fundamental difference between the mixed- 
spin systems with half-integer spins placed on sub-lattice B and respectively, 
their integer spins counterparts, consists in an appearance of the disordered 
phase in the former systems due to a presence of the 'non-magnetic' state 
S B = below D C /\J\ = —1.5. The latter systems, on the other hand, remain 
long-range ordered even below their last critical value D C /\J\ = —0.75, where 
their lowest but 'magnetic' spin state S B = 1/2 becomes stable. 

Now, we turn to a detailed investigation of finite-temperature phase diagrams 
as depicted in Fig. 2A (Fig. 2B) demonstrating the variations of critical tem- 
peratures with the uniaxial single-ion anisotropy for several values of half- 
integer (integer) spins constituting the sub-lattice B. In these figures, the 
critical temperatures of the mixed-spin Ising model on the 4-8 lattice (solid 
lines) are compared with critical lines of the mixed-spin Ising model on the 
honeycomb lattice (dotted lines) that represents another exactly soluble three- 
coordinated planar Ising lattice [4] . It can be clearly seen from Fig. 2 that both 
the mixed-spin Ising lattices exhibit very similar dependence on the single-ion 
anisotropy strength as the critical temperatures of the former lattice are merely 
slightly reduced with respect to those of the latter one. This result provides a 
further confirmation of the conjecture, originally supposed from a comparison 
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of critical points of the spin- 1/2 Ising model on the square and kagome lat- 
tices [11], that the critical points of Ising lattices with the same coordination 
number but different lattice topology being slightly lower for the lattices com- 
posed of more regular polygons. Finally, it is worthy to notice that the critical 
temperatures shown in Fig. 2A correctly reproduce an exact critical point of 
the isotropic spin-1/2 Ising model on the 4-8 lattice [12] by considering the 
limiting case D/\J\ — > — oo. 

For easy reference, let us introduce the classification scheme diagrammatically 
shown in Fig. 3 before proceeding to a discussion of the most interesting results 
obtained for temperature dependences of the resultant magnetization. It seems 
unavoidable that the Neel's nomenclature 1 referring to ferrimagnetic m(T) 
curves should be furnished by several less frequent but notable dependences. 
In addition to the familiar Q-, R-, P- and N-type curves, which are commonly 
found in the standard textbooks on magnetism [14], the L-, M-, S- and W-type 
curves are added here to this basic classification scheme. The conventional 
Q- and R-type dependences show a monotonic decrease of the magnetization 
with increasing the temperature, the former one displaying a steep decrease 
of the magnetization just in the vicinity of critical temperature, while the 
latter one already exhibits a relatively rapid magnetization decline within 
the range of intermediate temperatures before going abruptly to zero at the 
critical point. Another familiar P-type dependence shows the temperature- 
induced maximum as the temperature raises, whereas the N-type curve is being 
characterized by one compensation point at which resultant magnetization 
disappears due to the complete cancellation of the sub-lattice magnetization. 
The L-type curve [15] is very analogous to the P-type dependence, however, the 
resultant magnetization starts from zero in this particular case. Similarly, the 
M-type dependence also begins from zero but it shows two separate maxima 
before reaching the critical temperature. The stair-like S-shaped dependence 
initially exhibits a steep decrease of the magnetization that almost completely 
diminishes in the range of intermediate temperatures and finally, a second 
steep decrease follows near the critical temperature. Last, the striking W-type 
curve peculiarly exhibits two compensation points before the critical one. 

The temperature dependences of resultant magnetization normalized per one 
site of the mixed-spin 4-8 lattice (m = \m,A + ^b|/2) are depicted in Fig. 4. 
The total magnetization is plotted against the temperature for several values 
of the quantum spin number S and various single-ion anisotropies. At first, let 
us look more closely on the particular cases when the sub-lattice B consists of 



1 L. Neel [13] originally envisaged 6 distinct types of m(T) curves, but some of them 
having a non-zero initial slope in contradiction with the m(T) curves commonly 
found in the ferrimagnetic Ising systems. Therefore, the nomenclature used further 
is based on a set of the standard Q-, R-, P- and N-type curves [14] filled out by L-, 
M-, S- and W-type dependences. 
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the integer spin-S* atoms. As far as the mixed spin-1/2 and spin-1 system is 
considered, the standard Q-type curve observable for D > gradually changes 
into the R-shaped dependence by reversing the sign of single-ion anisotropy 
term. The shape of the R-type dependence can be attributed to more rapid 
thermal variation of the sub-lattice magnetization mg, which is thermally 
more easily disturbed than ttia under the condition D < 0. The situation 
becomes even a little bit more involved when assuming the mixed spin-1/2 
and spin-2 system. Before entering the disordered phase at D C /\J\ = —1.5, 
this system exhibits in the ground state another first-order phase transition 
between two spontaneously ordered phases having the spin-2 atoms in the 
states S B = 2 and 1, respectively. Consequently, slightly below the critical 
value D C /\J\ = —0.5 that determines the coexistence of both the ordered 
phases one also encounters the P-shaped dependence with the temperature- 
induced maximum arising from S B = 1 — > 2 thermal excitations (see the 
curve labeled D/\J\ = —0.6). By contrast, the opposite S B = 2 ^ 1 thermal 
excitations are responsible for the occurrence of the S-shaped dependence 
when selecting the single-ion anisotropy strength slightly above this critical 
value (see for instance D/\J\ = —0.45). 



Before concluding, the most remarkable features of the mixed-spin 4-8 lattice 
with the half-integer spin-5 1 atoms on the sub-lattice B will be summarized. 
It should be stressed that the systems also exhibit the Q-, R-, P- as well as S- 
shaped curves when varying the single-ion anisotropy strength (see the right 
panel of Fig. 4). The curvature of these dependences is apparently caused 
by the same reasons as discussed earlier for their integer-spin counterparts 
and hence, the discussion concerning with their origin being here omitted for 
brevity. It is worthy to recall, however, that the mixed-spin systems with half- 
integer spin-S* atoms possess the antiferromagnetically ordered ground state 
below their last critical value D C /\J\ = —0.75. On account of the preferable 
thermal excitations S B = 1/2 — > 3/2 and S B = 1/2 — > 5/2, respectively, the 
L-type curve develops for single-ion anisotropies close enough to this critical 
value (e.g. D/\J\ = —0.8 and —1.0). Notice that the L-shaped dependence 
should not be confused neither with the magnetization of systems exhibiting 
reentrant transitions (in our case, the total magnetization rises much more 
steadily near the zero temperature), nor with the m(T) curves of ordinary 
antiferrromagnets placed in a non-zero external magnetic field. Finally, it is 
noteworthy that the mixed-spin systems having the half-integer spin-S* atoms 
on the sub-lattice B indeed behave, due to the almost complete cancellation 
of antiferromagnetically arranged spins, as ordinary antiferromagnets when 
D/\J\ < —1.5. In this respect, the critical frontiers displayed in Fig. 2A should 
be better considered in this range of parameters as the critical boundaries 
determining disappearance of the spontaneous 'antiferromagnetic' order. 
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4 Concluding Remarks 

In the present article, the star-triangle mapping is being used for computing 
the exact phase diagrams and magnetization of the mixed spin- 1/2 and spin- S 
(S > 1) Ising model on the bathroom tile (4-8) lattice. The particular attention 
was focused on the influence of uniaxial single-ion anisotropy on the criticality 
and the shape of the magnetization vs. temperature dependence. Depending 
on the quantum spin number S and the single-ion anisotropy strength, the 
temperature variation of resultant magnetization has been found to be ei- 
ther of Q-, R-, P-, L- or S-type. The N- and W-type curves with one or two 
compensation temperatures, respectively, cannot appear in the system under 
investigation on behalf of the uniaxial single- ion anisotropy effect. 

It is worthwhile to mention that the relative ease of the calculation procedure 
presented here enables further remarkable extensions. Actually, the approach 
based on the star-triangle mapping transformation can appreciably be devel- 
oped to account also for: i) different nearest-neighbor interactions assigned 
to the bonds within the square plaquettes of the original 4-8 lattice and re- 
spectively, to the links between them; ii) second-neighbor interactions between 
the A-sites; iii) four-spin interaction between each B-site and its three nearest- 
neighboring A-sites; iv) the biaxial single-ion anisotropy acting on the 5-sites. 
Notice that the last modification necessitates following the extended version of 
mapping procedure worked out previously for the mixed spin-1/2 and spin-S' 
(S > 1) Ising model on the honeycomb lattice [16]. In the spirit of results re- 
ported on formerly (see e.g. Refs. [5] and [16]), it is quite reasonable to suppose 
that the extensions ii)-iv) may significantly modify the magnetic behavior of 
studied systems and leastwise the N- and W-type curves should be expected 
to emerge in these modified versions of mixed-spin Ising models. 
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Figure captions 

Fig. 1 Diagrammatic representation of the mapping scheme. The isotropic mixed- 
spin bathroom tile (4-8) lattice consisting of two inequivalent sub-lattices A 
and B constituted by the spin- 1/2 (empty circles) and spin-S* (solid circles) 
atoms, respectively, is mapped on the anisotropic spin- 1/2 Ising model on 
the Shastry-Sutherland (orthogonal-dimer) lattice with the intra- and inter- 
dimer interactions schematically illustrated as solid and dashed lines. 

Fig. 2 Phase diagrams of the mixed-spin Ising model on the bathroom tile (solid 
lines) and honeycomb (dotted lines) lattices. Figure on the right (left) illus- 
trates critical frontiers of the mixed-spin systems composed of the integer 
(half-integer) spin-S* atoms. The empty circles demarcate special critical 
points of the mixed-spin 4-8 lattice that correspond to the coexistence of 
two different magnetically ordered phases. 

Fig. 3 Diagrammatic representation of various types of the ferrimagnetic magne- 
tization vs. temperature dependences. Acronyms for particular cases are 
given by the circled letters. 

Fig. 4 The temperature variation of resultant magnetization of the mixed-spin 4-8 
lattice exposed for several values of quantum spin number S and diverse 
uniaxial single-ion anisotropy strengths. 
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5*= ±1/2 ms?=-S f - 




